Vortex dynamics in a three-state model 
under cyclic dominance 
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The evolution of domain structure is investigated in a two-dimensional voter model with three states under 
cyclic dominance. The study focus on the dynamics of vortices, defined by the points where the three states 
(domains) meet. We can distinguish vortices and antivortices which walk randomly and annihilate each 
other. The domain wall motion can create vortex-antivortex pairs at a rate which is increased by the spiral 
formation due to cyclic dominance. This mechanism is contrasted with a branching annihilating random walk 
(BARW) in a particle-antiparticle system with density-dependent pair creation rate. Numerical estimates for 
the critical indices of the vortex density (fi = 0.29(4)) and of its fluctuation (7 = 0.34(6)) improve an earlier 
Monte Carlo study [Tainaka and Y. Itoh, Europhys. Lett. 15, 399 (1991)] of the three-state cyclic model in 
two dimensions. 
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The self-organizing domain structures in the cyclic 
variants of the Lotka-Volterra model have been ex- 
tensively investigated because similar spatio-temporal os- 
cillations can appear in chemical reactions as well as in 
more complex ecological processes. These phenomena 
can be well studied within the formalism of voter mod- 
els H we will follow henceforth. Different versions of 
the three-state voter models under cyclic dominance on a 
square lattice were introduced by Itoh and Tainaka [^-^| . 
Bramson and Griffeath [Q studied the flux and fixation in 
cyclic particle systems. The pattern formation in cyclic 
cellular automata was investigated by Fisch Jy. Using 
a pair approximation, Frachcbourg and Krapivsky |9|] 
have shown that fixation occurs in a cyclic Lotka-Volterra 
model if the system is started from a random initial state 
and the number of states exceeds a critical value depen- 
dent on the dimension. According to this result, the sys- 
tem tends toward a self-organizing, inhomogeneous state 
if the number of states is less than 14 on a square lat- 
tice. Unfortunately, the theoretical understanding of the 
mechanism maintaining the inhomogeneous state is still 
incomplete. 

In this work we consider the features of a self- 
organizing domain structure in a two-dimensional system 
using the concept of vortices defined for three-state mod- 
els. Instead of studying the average sizes of the domains, 
Tainaka and Itoh H have determined the average density 
of vortices (c = (N v )/L 2 where N v denotes the number 
of vortices in a system with L x L lattice points). The 
vortices are points in these domain structures where the 
three different states (A, B and C) and the three types of 
domain walls meet. Evidently, the value of c _1 represents 
roughly the average domain area (size). The investiga- 
tion of vortex density was strongly motivated by the fact 



that its determination is much easier than the evaluation 
of the average domain size. 




FIG. 1. A typical domain structure for a three-state (color) 
model. The black bullets with white border and the white bul- 
lets with black border represent the vortices and antivortices. 

Figure 1 shows a (three-color) domain structure on the 
macroscopic scale. It is easy to recognize that two types 
of vortices may be distinguished as indicated by black and 
white bullets in the figure. We will call them vortex and 
antivortex depending on whether we find ABC or ACB 
order when going clockwise around the center. A simple 
rule may be deduced at first glance, namely the number of 
vortices is even around a closed domain. More precisely, 
the vortices and antivortices are alternately located along 
the closed boundary of a domain. This feature has serious 
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consequences when the motion and collision of vortices is 
considered. 

During the time evolution of a three-color domain 
structure the vortices move together with the bound- 
aries. In these processes the vortices can collide and 
annihilate each other. Figure || illustrates the typical 
elementary events (after the transient processes) whose 
combinations describe all the possible phenomena related 
to the creation, annihilation and collisions of the vortices. 
In the present approach the topological situations where 
four (or more) domain boundaries meet at a given point 
are considered as instantaneous events of collisions, pair 
annihilations or creations. Evidently, the illustrated ele- 
mentary events modify the connectivity among the vor- 
tices. This connectivity, however, can be modified by 
either fusion or fission of domains without any change in 
their constellation. 




FIG. 2. Shematic plots for characterizing the elements of 
vortex dynamics on a three-state map. The upper process (a) 
represents the annihilation of a vortex-antivortex pair walk- 
ing along their common boundary line as well as the reverse 
phenomenon corresponding to a spontaneous pair creation. 
(b) The "collision" of two vortices or antivortices can modify 
the type of domain separating them, (c) A vortex-antivortex 
pair can annihilate each other in a different way when crossing 
through the separating domain. 

A numeric analysis of the vortex density was performed 
by Tainaka and Itoh || in a two-dimensional voter model 
where the voters, located on a square lattice, should 
choose among three states (opinions): A, B and C . The 
system evolution is governed by a simple algorithm: a 
randomly chosen voter or one of its nearest neighbors 
can modify their opinion if those are different. If the 
chosen voter and its neighbor are in the states A and B, 
then the first voter adopts opinion B with probability 
P, otherwise its neighbor changes its state to A. This 
adoption rule is repeated cyclically for the B-C and C-A 



cases. Similar cyclic dominance characterizes the "paper, 
scissors, stone" games [ [Tcfl . 

The direction of dominance can be reversed by replac- 
ing 1 — P for P, therefore the analyses are restricted 
to P > 1/2. For P — 1/2 the present model is equiv- 
alent to a traditional voter model j|] which exhibits a 
(three-color) domain coarsening phenomenon if initially 
the voter states are random. In other words, the finite 
system evolves into one of the three homogeneous states 
while the vortex density goes to zero. 

For P > 1/2, however, a self-organizing domain struc- 
ture is maintained, and the vortex density tends to a 
stationary value, N v dependent on P. More precisely, 
Tainaka and Itoh have found a power law behavior, 
namely, 

c ~(p-p c r (i) 

where P c = 1/2 and f3 ~ 0.40 §. 

We have repeated these simulations using larger sys- 
tem size (400 x 400) and longer sampling times. In the 
vicinity of P c (P — P c < 0.01) the thermalization is 
chosen to be longer than 250,000 MCS (Monte Carlo 
Steps per particle). In the stationary state we have 
determined all the four-point configuration probabilities 
Q(n\, ri2, ri3, n^) [rii = A, B or C; (i = 1,2,3,4)] of a 
2x2 cluster. In our notation the (A, A, B, C) configura- 
tion refers to a vortex, the (A, A, C, B) to an antivortex, 
and the (A, B, C, A) configuration can be interpreted as 
a vortex-antivortex pair (before their annihilation or af- 
ter their birth). Monitoring the vortex density we were 
able to determine its average value and its fluctuation 
[x = L 2 ((N v / L 2 — c) 2 )] in the stationary states as a func- 
tion of p = P — P c . 
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FIG. 3. Log-log plot of vortex density vs P in the 
three-candidate voter model under cyclic dominance. The 
open diamonds represent MC data, the solid and dashed lines 
(resp. slopes 0.29 and 0.41) indicate the fitted power laws. 

The results of our simulations are summarized in two 
log-log plots (see Figs. || and ||). Figure || demonstrates 
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that we have reproduced Tainaka and Itoh's data H for 
p > 0.1. In this region the fitted power law (dashed line) 
is characterized by an exponent j3 = 0.41(3) in agreement 
with Tainaka and Itoh. For smaller p values, however, 
we have found a different exponent j3 = 0.29(4) (solid 
line). Similar crossover behavior can be observed for the 
fluctuation as shown in Fig. ^. The fluctuation remains 
approximately constant for p > 0.1, whereas it follows a 
power law (x ~ (P — Pi)~ 7 with 7 = 0.34(6)) for smaller 
p values. 




FIG. 4. Log-log plot of density fluctuation vs P — P c in 
the three-candidate voter model under cyclic dominance. MC 
results are represented by open squares, the fitted power law 
function is represented by a solid line for small P — P c values. 

The above mentioned configuration probabilities sat- 
isfy the reflection, rotation of 90° and cyclic symmetries. 
Consequently, all these quantities are describable, within 
a cluster approximation by introducing 6 independent pa- 
rameters which can be evaluated by solving a set of equa- 
tions of motion for the four-point configuration probabili- 
ties. Previously, this method has been proved to be a very 
efficient tool for the investigations of stochastic cellular 
automata p!l] |, evolutionary games p2[ , different lattice 
versions of the Lotka-Volterra models (at the level of pair 
approximation) |^,[|, and for a two-dimensional driven 
lattice gas maintaining a self-organizing domain struc- 
ture too p3| . In the present case, the calculations show 
a very weak P-dependence of the configuration probabil- 
ities including the vortex constellations mentioned above 
in contrast with the MC simulations. At the level 
of the pair approximation the configuration probabilities 
are independent of P ||, while the mean-field (1-point) 
approximation predicts homogeneous oscillatory behav- 
ior. At the five-point level, on the other hand, the prelim- 
inary calculations (integrating numerically the equations 
of motion) also indicate a very weak P-dependence. This 
puzzling failure of the dynamical cluster technique in- 
spired us to search for a mechanism observable at macro- 
scopic (or mesoscopic) scale. 



For P = 1/2, the domain coarsening is accompanied 
by the annihilation of pairs, and this process is not pre- 
vented by the weak spontaneous pair creation. Under 
cyclic dominance (P > 1/2), however, we have detected 
the appearence of a pair creation mechanism which is 
able to compensate for the previous annihilation process, 
yielding a finite vortex density. We have displayed in 
Fig. the evolution of a single vortex whose geometry 
allows us to recognize the essential processes. Here, in- 
stead of the usual periodic boundary conditions, we have 
assumed that a voter residing on the perifery and its "hy- 
pothetic" outer neighbor have always the same opinion. 
Along the boundaries one can observe a cyclic invasion 
whose average velocity component perpendicular to the 
border line is proportional to P — P c . This motion of 
the boundaries yields a spiral formation around the cen- 
ter. The average time evolution is decorated by noise as 
shown in a series of snapshots in Fig. || Due to the 
randomness the neighboring boundaries can contact and 
create vortex-antivortex pairs inside the spiral because it 
consists of narrow "arms" which is a favoured situation 
for the pair creation. The created pairs can be consid- 
ered as the offspings of the original vortex. Most of the 
pairs are annihilated within a short time but, sooner or 
later, a vortex-antivortex pair will eventually drift apart 
from each other. The corresponding vortices will then ex- 
pand and become capable to create further offsprings via 
the same spiral formation process. The self-organizing 
domain structure can thus be maintained by this mech- 
anism. 
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FIG. 5. Time evolution of a vortex initially having straight 
border lines for P = 1. The figures at the upper-left corners 
indicate the time measured in MCS units. 

Within the framework of the vortex language, the evo- 
lution of domain structure can be described by the two- 
dimensional motion of vortices allowing the annihilation 
and creation of pairs. Assuming that the motion of vor- 
tices is dominantly controlled by noise, the present prob- 
lem can be considered as a so-called branching annihi- 
lating random walk (BARW) with two types of particles 
created and annihilated in pairs. The branching process 
is evidently controlled by the value of P — P c , though the 
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mathematical relation between the branching rate and 
P — P c is not yet clarified. Furthermore, the branching 
rate is affected by the nearest neighbor distances because 
it limits the spiral formation. 

The traditional BARWs have been intensively studied 
in the last years (for recent reviews see the work by Cardy 
and Tauber jlj] and Marro and Dickman p5[ ) because 
they undergo a critical transition when varying the rate 
of branching. The corresponding critical behavior be- 
longs to the DP (directed percolation) universality class 
p6| involving the Reggeon field theory [|l7|], the surface 
reaction |l8| and Schlogl models and the extinction 
phenomena observed in spatial evolutionary games [ |r2| . 
According to the "DP conjecture" [^0) most of the one- 
component model with a single absorbing state belong to 
the DP universality class. Exceptions can appear when 
additional symmetries or conservation laws are intro- 
duced. Well known examples are those models in which 
the parity of particles is conserved during the elemen- 
tary processes |22|,^3| . The introduction of two (or more) 
types of particles has enlarged the number of possible 
universality classes |l4|,24 



The field theoretical results [[14 24 25j] indicate that 
mean-field approaches can give satisfactory description 
of the two-dimensional BARW models with two types of 
particles. This observation motivated us to try to de- 
scribe the branching annihilating random walks of vor- 
tices and antivortices through a mean-field equation with 
a density-dependent branching rate, namely 



d 

Ft' 



A(c)c 



(2) 



where c denotes the concentration of vortices and antivor- 
tices. The first term describes the annihilation process 
whose prefactor is eliminated by choosing a suitable time 
scale. For simplicity we suppose that the density depen- 
dence of the branching rate follows a power law with 
exponent v 



A(c) = pc 



(3) 



where p = A(P — P c ). Notice that this branching rate di- 
verges in the limit c — > if v > 0. In the stationary state 
one can easily determine the concentration as a function 
of p: 



(4) 



Comparing this formula with the above MC results one 
can conclude that the present mean-field description pre- 
dicts v ~ 1.5 if P — P c > 0.1 and v ~ 2.5 for the smaller 
value of P — P c . 

In order to check the role of fluctuations we have per- 
formed MC simulations on a particle-antiparticle BARW 
model. The system evolution is governed by nearest 
neighbor jumps, particle-antiparticle pair annihilations 
and creations as follows. A randomly chosen particle 
(or antiparticle) can create a particle-antiparticle pair lo- 
cated on two randomly chosen nearest neighbor sites with 



a probability Pb r , otherwise this particle jumps to one of 
the nearest neighbor sites. The processes which would 
result in two particles (or antiparticles) residing on the 
same point are blocked. A particle-antiparticle pair is an- 
nihilated if they would stay on the same point as a result 
of the mentioned elementary processes. The branching 
rate Pf, r is determined for a given particle as a product 
Pbr = PR1R2R3, where Ri (i?2, R3) denotes the distance 
between the chosen particles and its first (second, third) 
neighbor antiparticles. The choice of this branching rate 
is motivated by the topological fact that, in the origi- 
nal voter model, a vortex is connected directly to three 
antivortices by boundary lines (see Fig. |l|). During the 
simulations the value of Pb r can become larger than 1 
very rarely thus we did not need to reduce the time unit 
in which each particle has a chance to create offsprings or 
to jump. Initially particles and antiparticles with equal 
numbers are distributed randomly on a square lattice. 
The system size is varied from L = 100 to 400 when de- 
creasing p. During the simulations we have recorded the 
number of particles and determined the average value of 
the concentration c and its fluctuation \ in the station- 
ary state as defined above. From these MC results (see 
Fig. ^) we could confirm that the concentration follows a 
power law (solid line) with an exponent — 0.42(3). It is 
emphasized that this value of f3 agrees very well with the 
prediction of the above mean-field formalism (Pmf = 0.4 
for v = 1.5). This result implies that any value of the ex- 
ponent P is reproducible with the parameter adjustment 
of a more sophisticated BARW model. 




FIG. 6. MC results for the particle-antiparticle concentra- 
tion (open diamonds) and its fluctuation (open squares) as a 
function of p in the BARW model described in the text. The 
solid line (slope 0.42) indicates a fitted power law. 

Considering the fluctuations obtained by simulations, 
a striking difference is found between the present BARW 
model and the three-candidate voter model. The BARW 
simulations (squares in Fig. ||) indicate that x{p) ^ c (p)t 
which seems to be a typical behavior for the BARW mod- 
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els p5| , p6| . For the voter model, on the contrary, \ di- 
verges (see Fig. ^) for small values of the control para- 
mater. We can identify two possible sources for this dis- 
crepancy. First: during the diffusive motion the particle- 
antiparticlc annihilation results in an aggregation of the 
two species [^7t , whereas in the vortex-antivortex system 
such a process is strongly limited by the mentioned topo- 
logical features. Second: the three- fold degeneracy of the 
absorbing states of the three-candidate voter model has 
no counterpart in the traditional BARW models. Notice 
that, in the suggested BARW models, the evolution into 
an absorbing state is prevented by the divergency of the 
branching rate in the limit c — > for sufficiently large 
system size. Further systematic research is required to 
clarify the effect of these phenomena. 

To summarize, in the present paper we have improved 
the accuracy of the numerical analysis of the critical 
transition appearing in the vortex density for the three- 
candidate voter model when varying the magnitude of 
cyclic dominance. Recognizing that the dynamics of the 
vortices and antivortices is similar to a BARW model 
with a density-dependent particle-antiparticle pair cre- 
ation (branching) rate, we have contrasted these two sys- 
tems. According to our comparison, we can state that the 
power law behavior of the vortex density is reproducible 
with a suitable choice of the pair creation mechanism. 
The same is not true for the behavior of fluctuations, 
which seems to be quite different in the two models. This 
discrepancy is a motivation to seek further extensions of 
BARW models, since this approach seems to be very use- 
ful in the investigation of the self-organizing, three-color 
domain structures. 
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